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Abstract 

We discuss the interplay between lepton asymmetry L and v oscillations in the 
early Universe. Neutrino oscillations may suppress or enhance previously existing 
L. On the other hand L is capable to suppress or enhance neutrino oscillations. The 
mechanism of L enhancement in MSW resonant u oscillations in the early Universe 
is numerically analyzed. L cosmological effects through u oscillations are discussed. 
We discuss how L may change the cosmological BBN constraints on neutrino and 
show that BBN model with z/e -f-)- Us oscillations is extremely sensitive to L - it 
allows to obtain the most stringent constraints on L value. We discuss also the 
cosmological role of active-sterile z/ mixing and L in connection with the indications 
about additional relativistic density in the early Universe, pointed out by BBN, 
CMB and LSS data and the analysis of global u data. 
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Introduction 

Lepton asymmetry (L) of the Universe is not measured yet and may be orders of magnitude 
bigger than the baryon asymmetry /3 ~ 6.10~^°, which was measured with great precision 
from Cosmic Microwave Background (CMB) and Big Bang Nucleosynthesis (BBN) data. 
L is usually defined as L = (A*"; — Ni)/N^, where Ni is the number density of leptons, 
Ni of antileptons, while is the number density of photons. Considerable L might 
be contained only in u sector. Thus, the detection of the Cosmic Neutrino Background 
would provide L direct measurement. Till then L is measured indirectly by its influence 
on observable relics of the Universe. The abundances of the primordially produced light 
elements during BBN provide such a sensitive test of L. 
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Main cosmological effects of L on BBN are its dynamical and Icinetic effect. The dy- 
namical effect of L consists in tlie increase of tlie radiation energy density of tlie Universe 
due to non-zero L: 5pr = [1 + 7 /8{4/liy/^6Neff]p^, wliere AN.ff = 15/7[(^/7r)^+2(e/7r)2], 
wliere ^ = n/T is the v degeneracy parameter. Tliis effect leads to faster Universe expan- 
sion B. = (S/STrGpy^"^, which leads to earlier freezing of nucleons and hence influences 
BBN produced elements. The kinetic effect is noticeable for big enough \L\ > 0.01 in 
the sector and is due to different number densities of Ue and u^, which participates in 
neutron- proton transfers: + n p + e" , 17^+ p n + e'^ in the pre-BBN epoch and, 
correspondingly, BBN yields (see ref. [1] and refs there in). 

Besides these well known effects L may influence BBN through its interplay with 
neutrino oscillations. Unlike previous effects, this L effect on BBN, called further indirect 
kinetic effect, may be noticeable for much smaller L values, L « 0.01 [2]. Hence BBN 
with oscillations provides the possibility to measure and/or constrain tiny L with values 
close to /3. 

In the next section the interplay between L and u oscillations is described. In the third 
section the indirect kinetic effect of L on BBN is studied and the possibihty to constrain 
L via BBN with u oscillations is discussed. L role as a solution of the problem of the 
excess radiation density in the Universe is provided in the fourth section. 

1 Interplay between L and neutrino oscillations in 
the early Universe 

Small L << 0.01 has negligible dynamic and direct kinetic effect. Nevertheless, due 
to its interplay with u oscillations such L is capable of changing u number densities, u 
distribution and spectrum distortion and changing u oscillations pattern (suppressing or 
enhancing oscillations), by which it influences nucleons kinetics and flnally BBN produc- 
tion of hght elements [2]. This effect was proven to persist down to L ~ 10~^ [2, 3, 4]. 

We studied numerically the interplay between tiny L, 10~^° < L < 10 and electron- 
sterile Ue ■<-)■ oscillations, effective after active u decoupling 5m^sin^2^ < 10~^ eV^. 

It is known that active-sterile oscillations may change neutrino- antineutrino asymme- 
try of the medium, suppress or enhance it [2, 3, 4, 5, 6, 8, 9, 10, 11, 12, 13, 14, 15, 16]. 
^ On the other hand L influences v propagation. Qualitatively this influence may be 
described as follows: The average potentials Vf for v depend on the particle asymmetries 

^ There are other cosmological effects of active-sterile neutrino oscillations, like excitation of additional 
light particles into equilibrium [5, 6] and distortion u energy distribution [7, 8]. 
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of different constituents of the medium and they differ for different neutrino types due to 
different interactions with the particles of the plasma [17]: 

Vf = V2GfQN^/M^ ± LN^ (1) 

where / = e,iJ,,T, "minus" corresponds to u and "plus" to u, Q ^ —ET, L ~ —La, 
L" is given through the fermion asymmetries of the plasma (in the discussed case L ~ 

In the adiabatic case the effect of the medium can be hidden in the oscillation param- 
eters Sm^ and by introducing matter oscillation parameters: 

sin^^m = sin^ ^/[siii^ ^ + 2Gf{{Q/M^ ^ L)N^/5m^ - cos2'df]. (2) 

L, as a characteristic of the medium, may suppress oscillations by decreasing their ampli- 
tude, or enhance oscillation transfer in case resonant condition between the parameters 
of the medium and the oscillation parameters holds: 

2Gf{Q/M^ ^L)N^ = cos2i? X Sm"^ (3) 

In the early Universe at high temperature Q/M^ > L resonant oscillations both for u and 
P are possible if Srn^ < 0. With the cooling of the Universe when L begins to dominate, 
Q/M^ < L, resonant transfer for antineutrinos in case Sm^ < 0, or for neutrinos if 
5m^ > is possible. ^ Thus neutrino propagation and resonance in the neutrino sector 
differs from that of antineutrino for non-zero L. Due to L influence of the u propagation 
L may change rii,, its spectrum distribution and oscillation pattern. 

This simplified description of the medium influence is applicable in the equilibrium 
situation, when working in terms of average momentum and particle densities is reason- 
able. In the nonequilibrium situation, when spectrum distribution of p is considerable, 
as is the case of late electron-sterile oscillations discussed here, it strongly effects both u 
propagation and L evolution. Hence, for the correct description of the neutrino - asym- 
metry interplay it is essential to provide an accurate account of the neutrino spectrum 
distortion due to L and oscillations. ^ 

^The resonant condition (3) for Q = was first studied and is known as Mikheev-Smirnov-Wolfenstein 
effect [18]. 

^For example, in the nonequilibrium situation, when spectrum distribution of u was properly described, 
simultaneous resonance transfer was found possible also in the Sm^ < case due to the "resonant wave" 
passing through the neutrino distribution [2]. 
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1.1 Exact description of the propagation of neutrinos and L 
evolution 

The equations governing v evolution, in terms of neutrino density matrix in momentum 
space, are given below. They account simultaneously for neutrino-L interplay. Universe 
expansion, v oscillations and v forward scattering and describe precisely v energy distri- 
bution [8]. Our numerical analysis of L and v propagation was based on these equations: 

dp{t)/dt^Hp, {dp{t)/dp,) + 
+i [Ho: P{t)] + iV2GF {C - Q/M^) [a, p{t)] + O (G^) 
dp{t)/dt = Hp, {dp{t)/dp,) + 
+i [Ho, p{t)] + i^GF {-C - Q/M^) [a, p{t)\ + O (gI) . 

Lu, ~ / d^p(pLL - Pll)/N^ 



where «jj = U*^Uje, Vi = UiiVi{l = e,s). Ho is the free z/ Hamiltonian. Q arises as an 
W/Z propagator cffect,g ^ T. 2L,^ + L,^ + L^^, L^,^ ~ (A^^,^ - Np^f)/N^. At 
decoupling of z/g Ug was assumed empty. 

Due to L the equations arc coupled integro-diffcrcntial. L leads to different evolution 
of u from z/ due to the different sign of L in the equations. Numerical analysis of the 
evolution of u ensembles, evolution of L, and also the evolution of nucleons for the entire 
range of oscillation parameters and for the temperature range [2 McV, 0.3 MeV] and for 
10~^° < L < 0.01 was provided in non-resonant u oscillations case. In case of resonant 
oscillations L initial value was taken to be Lj ~ /3. We have described precisely u momenta 
distribution: 5000 bins were used in the non-resonant oscillations case, and up to 10 000 
in the resonant case. 

In case of nonresonant oscillations and relic L the following relations describe with 
good accuracy the exact behavior of L and L-oscillations interplay: L > 10~^ en- 
hances oscillations, L > 0.1(5m^/eV^)^/^ suppresses oscillations, and asymmetries L > 
{Sm'^/eV^y^^ inhibit oscillations. For illustration of the exact dependence see the figures 
in ref. [4]. 

In case of resonant Ve ■'^ oscillations the evolution of L has a rapid oscillatory 
behavior. The region of parameter space for which a generation of L is possible was found 

*The case of non-zero population of was considered in refs. [15, 19, 20]. 
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|5m^| sin^ 26 < 10^^'^ eV^. ^ A maximum possible growth of L by 5 orders of magnitude 
was determined. 

L role in BBN with v oscillations was numerically studied as well. The change in 
BBN constraints on oscillation parameters due to L and BBN constraints on L in case of 
Ve ^ are presented in the next section. 

2 BBN with active-sterile neutrino oscillations and 
lepton asymmetry 

Big Bang Nucleosynthesis is theoretically well established due to the precise data on 
nuclear processes rates, known from laboratory experiments at the energies relevant to 
BBN epoch, and the precise data on the light elements abundances D, He and Li. Besides, 
the baryon-to-photon ratio, which is the only parameter in the standard BBN, has been 
independently measured with good accuracy by CMB precision data. Therefore, BBN 
is used as the most early and precision probe for physical conditions in early Universe 
and, hence, presents the best astrophysical and cosmological probe for new physics and 
microphysics, relevant at BBN energies. 

Being the best speedometer at radiation dominated stage BBN was used to probe u 
properties, the number of light species, u oscillations, distortion in u distribution, etc. 
Due to its sensitivity to the expansion rate and to the nuclcons kinetics in the prc-BBN 
epoch it has been shown that BBN presents also the most exact leptometer (see ref. [4] 
and references there in). Thus BBN studies allow to put stringent limits on oscillation 
parameters [5, 7, 6, 8, 21, 12, 13, 15, 19, 20, 23] and on L in the presence of electron-sterile 
oscillations [2, 3, 4]. 

2.1 BBN constraints on L 

Eventual u degeneracies of flavor neutrinos will equilibrate before BBN due to flavor u 
oscillations, having in mind recently measured value ^13. Then BBN constraints on L read 
[24, 25, 26, 27, 28, 29]: 

\L\ < 0.1. (4) 

These constraints allow the possibihty of L orders of magnitude larger than the mea- 
sured /3. However, for such small L the dynamical cosmological effect of L isnegligible. 

^The instability region is slightly more stringent than the existing in literature for other oscillation 
models ref. [22]. 
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Hence, that small L in the flavor v sector cannot mimic the extra relativistic degrees of 
freedom during BBN, which seem to be required by recent analysis of cosmological and 
neutrino oscillations data (to be discussed in more detail in the fourth section). 

2.2 L and BBN with late electron sterile v oscillations 

In case oi L /3 stringent BBN constraints on active-sterile oscillation parameters ex- 
ist. The presence of L > (3 may change considerably BBN constraints on oscillations 
parameters due to L indirect kinetic effect. 

In BBN with electron-sterile v oscillations Yp decreases at small mixing parameters 
values due to L growth caused by resonant u transfer, thus at these mixing angles oscil- 
lations generated L relaxes the BBN constraints on oscillations [8, 12, 4]. 

In case of relic L > /3 the presence of such L relaxes the BBN bounds on neutrino 
oscillation parameters in case of electron-sterile non-resonant neutrino oscillations at max- 
imal mixing and strengthens them at small mixing angles. In this case depending on the 
asymmetry-oscillations interplay the asymmetry may enhance, suppress or stop u oscilla- 
tions, reflecting correspondingly to strengthening, relaxation or elimination of the BBN 
constraints. In the last case the approximate BBN constraint reads: 

W(ey2) < L=^/2. (5) 

Vice versa, this constraint can be considered also as a cosmological constraint on L from 
BBN with electron-sterile neutrino oscillations. Thus, L generated in the electron-sterile 
sector may (partially) suppress the oscillations in other sectors, hindering equalization of 
chemical potentials, thus relaxing the stringent BBN bound. The presence of L (no matter 
how generated), capable to suppress oscillations, may lead to only partial population of 
the sterile neutrino. Vise versa L values capable to enhance oscillations will bring faster 
the sterile state into equilibrium. 

The model of BBN with ^ i^s, cfTcctivc after electron u decoupling presents the 
possibility to feel extremely small asymmetries L > 10~^ due to the indirect kinetic effect 
of L (see ref. [4]). This sensitivity allows to derive stronger constraints on L. Namely, 
equation (5) can be used to put the following BBN constraint on L: Having the indications 
for active-sterile oscillations with Sm"^ ~ 10~^ [30] and replacing this value in eq.(5), a 
much stronger upper limit on L (than the one presented in eq.(4)) follows L < 10~^'^. 
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3 L and excess radiation density in the Universe 



In recent years an increasing number of cosmological indications suggesting excess rela- 
tivistic density, corresponding to different epochs, appeared [31, 32, 33, 34, 35, 36, 37]. ^ 
Besides, v oscillations data require 1 or 2 additional sub-eV sterile neutrino, participating 
into oscillations with flavor neutrinos with higher mass differences values, than the ones 
required by solar and atmospheric u oscillations experiments [40, 42, 43, 44, 45]. It is 
interesting if cosmology allows 2 light additional sterile neutrinos and if they can explain 
the excess relativistic density. 

Additional light sterile neutrinos with the mixing and mass differences estimated by u 
oscillations data with mass differences in the eV range will be broTight into equilibrium in 
the early Universe. BBN favors the presence of one such but He and D data excludes 
2 fully thermalized Ug. Besides, neutrinos in sub-eV range produce too much hot dark 
matter [46, 47]. Thus, two additional i/g are in tension both with BBN and with LSS 
requirements [48, 49, 36]. 

L, namely its dynamical and direct kinetic effects, has been considered as an expla- 
nation of the excess radiation. It was shown that excess radiation cannot be explained 
by degenerate BBN [29]. However, the presence of L may be the solution in case its 
value is enough to suppress activc-stcrilc oscillations so that Ug are not fully thermal- 
ized [6, 9, 2, 41]. Our estimation of the value of L necessary to suppress oscillations and 
achieve the suppression of Us production is L > 0.08. This is higher than the values 
discussed by refs.[41, 51], that found \L\ > 10^^. The difference might be due to different 
approximations used or to the fact that previous studies do not account precisely for the 
1/ energy distribution. 

Thus, in modifled BBN with u oscillations and high enough L the models with ad- 
ditional hght sterile neutrinos may be allowed. To obtain the exact L value a precise 
numerical analysis, solving the exact kinetic equations, including all u species and ac- 
counting for all L effects, discussed above, should be provided. 

Hence, the excess relativistic density might point to additional sterile neutrinos and 
the presence of L. However, there exist other possibilities as well, namely: u active-sterile 
late oscillations leading to the overproduction of He-4 and thus imitating extra radiation, 
MeV decaying particles during BBN [53], or other modiflcations of the standard cosmo- 
logical model. Future experimental and observational data will choose among different 
possibilities. In particular, it is expected that Planck data will be able to check with 

^See, however, the recent results [38], which are in agreement with the canonical value of Neff and 
also the discussion in ref.[39]. 
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higher sensitivity the status of extra radiation. 

4 Conclusions 

There exists interesting interplay between small L, relic or produced by active-sterile 
V oscillations, and late active-sterile v oscillations. A detail numerical analysis of this 
interplay between lepton asymmetry L « 0.01, and v oscillations was carried out. The 
evolution of v and L was studied using exact kinetic equations for v density matrix in 
momentum space and describing the v energy distribution with very high accuracy - up 
to 10 000 bins. 

A considerable enhancement of L - by 5 orders of magnitude - was found possible in 
non-equilibrium electron-sterile resonant neutrino oscillations, effective after the decou- 
pling of active neutrino. The region in the oscillation parameter space of considerable L 
growth was determined. 

The parameter ranges for which relic L is able to enhance, suppress or inhibit non- 
resonant electron-sterile oscillations was determined. 

Cosmological influence of small L, which do not have direct effect on nucleons kinetics 
during BBN, was discussed. Such small asymmetries are invisible by CMB at present, 
but may be felt by BBN. It was shown that L as small as 10~^ may be felt by BBN via 
oscillations. Also BBN constraints on v oscillations parameters depend nontrivially on L. 

L generated by oscillations at small mixing angles suppresses u oscillations and reduces 
the overproduction of Yp and relaxes BBN constraints. 

Relic L present during BBN, depending on its value, can strengthen, relax or wave 
out BBN constraints on oscillations. The value of L, capable to hinder i/ oscillations is 
determined for different sets of oscillations parameters of the model, a good approximation 
to the exact value is L > [bvp? j [eV'^yf'l^ . On the other hand this can be considered as the 
most stringent cosmological constraint on L based on BBN with electron-sterile neutrino 
oscillations. 

Due to its interplay with neutrino oscillations L may play important role for resolving 
the problem of additional radiation density suggested by cosmological observational data 
and by v oscillations data. In principal BBN with L may allow 3+1 and 3+2 oscillations 
models, in case L value is high enough to suppress active-sterile oscillations, thus providing 
incomplete thermalization of and relaxing BBN constraints on additional light sterile 
species. 
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